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^ ' Abstract 
<T 

For all non-symmetric discrete relativistic Toda type equations we establish a 
ON . relation to 3D consistent systems of quad-equations. Unlike the more simple and 

better understood symmetric case, here the three coordinate planes of Z 3 carry 
different equations. Our construction allows for an algorithmic derivation of the 
zero curvature representations and yields analogous results also for the continuous 
time case. 
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1 Introduction 



> 



c3 



This paper is devoted to one aspect of the general topic of discrete integrable systems. 
In the recent years the viewpoint that discrete integrable systems are in a sense more 



fundamental than continuous ones becomes gradually more and more accepted in the 
soliton theory community. The last developments led to the understanding that the very 
definition of integrability becomes more transparent and natural on the discrete level, 
see pj~]. One can consider certain types of discrete equations, like the so called quad- 
equations, as a sort of "elementary particles" of the solitonic world, with the enormous 
richness of this world resulting from various combinations and limiting procedures of 
these elementary objects. Here, we discover such "elementary particles" underlying the 
so called equations of the relativistic Toda type. This "microstructure" become visible 
5_i ■ only after the discretization procedure, and even then some hidden degrees of freedom 



remain to be uncovered in order for the whole simplicity to become apparent. 

We review the relativistic Toda type equations and their integrable discretizations in 
Sections[2]and[3l respectively. After that, general Toda equations on graphs are discussed 
in Section HI Fundamental notion of quad-graphs and quad-equations (which belong, 
in our view, to the "elementary particles" mentioned above) are discussed in Section 
[5l while the following Section [6] is devoted to the explanation of how the Toda type 
equations on graphs can be reduced to the systems of quad-equations. In Section [7] some 
combinatorial aspects of the regular triangular lattice are discussed, since this lattice 
underlies the discrete relativistic Toda type equations. Finally, Sections IHHTUI contain 
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the main new results of this paper. Namely, in Section [8] we present the systems of quad- 
equations which yield the non-symmetric discrete relativistic Toda equations. In Section 
[9] we show that all this systems have their common origin in just one master system of 
this type. In SectionQIJlwe demonstrate how to derive the zero curvature representations 
for relativistic Toda type equations (discrete and continuous) in an algorithmic manner. 
A brief outlook is formulated in the concluding Section [TT1 



2 Lattice equations of the relativistic Toda type 

The term "equations of the relativistic Toda type" is used to denote integrable lattice 
equations of the general form 



Xk = r(x k ) [x k+1 f(x k+1 - x k ) - x k _if(x k - x k -i) 

+g(x k +i - x k ) - g(x k - x k -i) )■ (1) 



The relativistic Toda lattice proper was invented by S. Ruijsenaars |14| . It is described 
by Newtonian equations of motion 

e Xk+l~Xk 

x k = (1 + ax k+ x)(l + ax k )- 



1 + a 2 e Xk + 1 ~ Xk 

gX k x k — 1 

-(1 + ax k )(l + aifc-i) — . (2) 

1 + a 2 e Xk ~ Xk - 1 

Here a is a small parameter whose physical meaning is the inverse speed of light. In the 
non-relativistic limit a — > system ([2]) turns into the usual Toda lattice. 

It took about a decade for further integrable equations of the relativistic Toda type 
to be discovered. In (16] the following ones were found: two systems which again can 
be considered as a-perturbations of the usual Toda lattice: 

x k = (l + ax k+1 )e Xk + 1 ~ Xk - (l + ax k -i)e Xk ~ Xk - 1 

_ a 2 e 2(x k+1 -x k ) + a 2 e 2{x k -x k ^) (3) 

and 

x k = (1 - a± k f ((1 - axk+i) e Xk +^ ~ Xk - (1 - ax k ^) e Xk ~ Xk -^, (4) 

and two systems which can be considered as a-perturbations of the so called modified 
Toda lattice: 

x k = x k (e Xk ^- Xk - e xk-x k -i 



I e x k +i-x k e x k -Xk-i \ 

+a x k+ iik Xkik-i (5) 

\ l + ae Xk + 1 ~ Xk l + ae Xk ~ Xk ~ l J 
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and 



/ e x k+1 -x k e Xk~X k -l \ 

x k = x k (l - ax k ) (1 - ax k+1 ) — - (1 - ax k ^) . (6) 

\ 1 + ae Xk+i-x k i + ae x k -x k -ij 

For all these systems a Lagrangian and a Hamiltonian formulations were given, and their 
complete integrability was demonstrated by presenting the full set of integrals of motion 
and a local zero curvature representations of the type 

L k = M k+1 L k - L k M k (7) 

in terms of 2 x 2 matrices (locality means that the matrix L k depends only on x k 
and the corresponding canonically conjugate momentum p k , and not on phase variables 
from other lattice sites). Two further systems of the relativistic Toda type with rational 
interactions (as opposed to exponential interactions in the previous ones) appeared in 

P2I= 

.. . , „ \ - ax k+ ix k ax k x k -i 

x k = x k (x k+ i - 2x k + x k -i) + — — - — — (8) 

l + a(x k+1 -x k ) 1 + a(x k - x k -i) 

and 

. /-. , 2 - \ f x k+i ~ x k - ax k+ i x k - x k -i - . . 

x k = x k (l + a x k ) r — . (9) 

V 1 + a(x k+ i - x k ) 1 + a(x k - Xk-i) J 

Both these systems are a-perturbations of the so-called dual Toda lattice. Also for 
these systems the Lagrangian and the Hamiltonian formulations were given in [T7], as 
well as a demonstration of their complete integrability. However, local zero curvature 
representations were not given at that point. 

Along another line of research, a complete classification of "integrable" systems of 
the type ([T]) was achieved in [6]. The notion of "integrability" used in this paper is 
a clever and unexpected device, allowing to carry out a complete classification, but it 
has, a priori, nothing to do with the usual Liouville- Arnold integrability. Namely, they 
noticed that the above systems are always Lagrangian, and required that their form be 
retained under a sort of Legendre transformation. This allowed them to find all the 
Newtonian equations mentioned above, as well a new series of systems, including the 
most general one, 

_l/. 2 _ 2^ ( sinh 2 ( x k+i ~ %k) ~ V 1 sinh(2^q) x k+1 

Xk ~ n\ X k u ) \ 2 . 2 

z \ sinh (x k+ i — x k ) — sinn [vol) 

sinh 2{x k — x k -\) — v~ x sinh(2i/a) x k -i \ 
sinh 2 (xfc — x k -\) — sinh 2 (z^a) J 

its limiting case (first rescale x k \— > vx k , v \— > 7^, and then send v — > 0): 



(10) 



/ -2 2\ ( x k+l ~ x k ~ ax k+l x k ~ x k-l ~ & x k-l \ . , , , 

Xk = -{X k -j) [- ^ T^2~l ^2 2~2 ' t 11 ) 

V {x k+ i - x k y - 'y a (x k - x k -iY - Ya z 



3 



as well as the particular cases v = 0, resp. 7 = 0, of the latter two systems: 
x k = -x\ coth(x fc+ i - x k ) - coth(x fc - x k -i) 



sinh 2 (x fc+ i - Xk) sinh 2 (a; fc - x k -i 
and 

■ ±k = _ ± 2 ( 1 <*Sfc+i + _axk-i_\ (13) 

k \x k+1 - x k x k -x k ^i (x k+1 -x k ) 2 [x k -x k _ x ) 2 ) 

Complete integrability of these systems in the usual sense, along with 2x2 local zero 
curvature representations for the last two ones, was demonstrated in the monograph 

EE]. 



3 Time discretization of the relativistic Toda type equa- 
tions 

Integrable discretizations of the relativistic Toda type equations all have the following 
general shape of discrete time Newtonian equations of motion: 

F{x k - x k ) - F(x k - x k ) = 

G{x k+ i - x k ) - G(x k - x k -i) + H(x k+1 - x k ) - H(x k - x fc _i). (14) 

Here and below we use the following abbreviations for functions of the discrete time KL: 

Xk = x k (t), x k = x k (t + h), x k = x k {t-h). 

The integrability preserving time discretization of the Ruijsenaars' relativistic Toda 
lattice was performed in [15j . where the following equations were derived: 

1 + ahrH e^k-Xk -i\ fi + a 2 e x k +i-Xk\ h + a ( a _ h ^ e x k -x k A 

) Y = 7 V7 T ' (15) 

1 + ah -l l e Xk~X k _ 1 \ M _|_ a 2 e X k -X k ^ \( 1 + a ( a _ h ) e Xk+l ~ X k \ 

It was shown in [15] that, upon a natural Lagrangian (or discrete Hamiltonian) re- 
formulation in terms of the canonically conjugate variables x k ,p k , discrete time equations 
(|15p share integrals of motion with the continuous time equations ([2]) , and therefore have 
the same integrability properties (belong to the same integrable hierarchy). This remains 
true also for all integrable discretizations in this section. 

The Newtonian systems ([3|) and Q were discretized in [16] in an "additive" manner 

as 

x k -x k _x k -x k = hae xk+i-x k _ haeXk - Xk _ 1 

h{a-h)e~ k+1 ~ Xk | h{a-h)e x k-xk-i 
1 - hae~ k+1 ~ Xk 1 - hae x k~Xk-i 



4 



and 

1 1 
l-ah- x ( e^k-Xk-i ] 1 - a h~ l ( e Xk ~~ k - 1 



/ i u\ x k+l ~ x k i , r\ x k~ x k-l 2 ^fc+l - ^fc . 2 x k~ x k-l /-, ~\ 

= a(a + nje — a(a + fije — a e + a e , (17) 

while the discretizations of the systems ([5]) and ([6]) given there were "multiplicative" : 
x k -x k _ ^ ( 1 + he x k+1 ~ x k^( 1 + aeXk -x k _^ ( 1 + ae x k+1 -x k 

(18) 

x k~Sk _ 1 ) l 1 + ae x k+i -Xk)l 1 + he x k -x k ^ + ae x k -x k ^ 



and 



> x k- x k-i) l -ah 



l / x k~ x k 



i Xk Sk_ 1 \ i_ a h-i\^ e x k -x k _i 

1 + ae Xk ~ Xk - 1 ) ( 1 + (a + h) e~ k+1 ~ Xk 



1 + ae x k+i~ x kj f 1 + (a + h) e Xk ~ x *- 

Discretizations of the rational systems © and © appeared in |17j : 

(l + h{x k+1 - x k )j (l + a(a; fc - x k -i)j (l + a{x k+1 - x k )j 



x k x k 

x k x k 

and 



1 + a(x k +i - x k ) ) ( 1 + h(x k - Xk-i) ) (l + a(x k -x k - 



-l 



(x k - Xk) (l + a(a + h)h~ 1 (x k - x k )j 
(x k - x k ) (l + a(a + tyh^ixk ~ x k)^j 

1 + a(x k -Xk-i)) [l + {a + h)(x k+1 -x k ] 



1 + a(x k +i -x k )) (1 + (a + h)(x k - x k -i) 



(19) 



(20) 



(21) 



A discrete version of the device from [6] was developed in [T] . It was used to classify 
"integrable" systems of the general type (I14p . i.e., those retaining their form under a 
sort of a discrete Legendre transformation. The resulting list consisted essentially of 
the systems quoted in above in this section, as well as of discretizations of the systems 
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(fID ]H [T5 ]) . A discretization of (OH reads 

sinh(xfc — Xk + vh) sinh(xfc — x& — vh) 
sinh(xj! c — Xk — vh) sinh(xfc — x& + vh) 

sinh(xfc + i — Xk + va) sinh(x£. — x&_i — ^a) 

sinh(xfc + i — Xk — va) sinh(xfc — Xk~\ + va) 

sinh( Xk+\ — Xk — va + vh) sinh(xfc — + va — vh) 

x . (22) 

sinh(xfc_|_i — Xk + va — vh) sinh(xfc — Xfc_i — va + vh) 

Its rational version, which is a discretization of reads 

(x k - x k + 7/1) (x k - x k - jh) (xk+i -x k + 7a) (x k - Xk-i - 7a) 
__ . = . x 

(x k -x k - jh) (x k - x k + 7/1) {x k+ i -x k - 7«) (^fc - Xk-l + 7«) 

(zfc+i - x k - 7« + 7^) 0^ - ^fc-i + 7« - lh) 

x . (23) 

(z fe+ i - x fc + 7a - 7/1) (x fc - x fc _i - 7a + 7/1) 

Finally, the v — > 0, resp. the 7^0 limits of the latter two equations lead to additive 
ones: 

h coth(xfc — Xk) — h coth(xfc — Xk) = a coth(xfc + i — Xk) — a coth(x/j — Xk-\) — 

— (a — h) coth(xfc+i — Xk) + (a — h) coth(x/c — Xk-i) (24) 

and 

h 

•^k Xk 

which discretize (|12|) and Q13p. respectively. The integrability of the difference equations 
(|22p ^ (|25p in the usual sense was dealt with in [18], where it was shown that they share 
the integrals of motion and the matrices Lk from the zero curvature representations with 
their continuous time counterparts. 

4 Discrete Toda type equations on graphs 

An important observation made in [2] was that the natural combinatorial structure un- 
derlying the discrete relativistic Toda type lattices (fTi|) is actually the regular triangular 
lattice in the plane (rather than the standard square lattice Z, 2 ). Namely, each equation 
of the system (|14p relates seven fields assigned to the star of a vertex of the regular 
triangular lattice, each one of the functions F, G, H being associated to edges of one of 
the tree directions, see Figure CD 

Note that while in the systems (|15p - (|2ip (called hereafter non-symmetric discrete 
relativistic Toda type equations) the tree functions F, G, H are, generally speaking, dif- 
ferent, this is no more the case for systems (|22p - (|25p . In the latter systems (called 
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a — h 



+ 



a — h 



Xk Xk Xk-l-l Xk Xk Xk — l Xk-i-\ Xk Xk Xk — l 



(25) 



Xk-1 x k 




•£k Sk+l 



Figure 1: Regular triangular lattice underlying discrete relativistic Toda type systems 



hereafter symmetric discrete relativistic Toda type equations) all three functions essen- 
tially coincide, differing only by the values of the built-in parameters. 

This allows, at least for symmetric systems, the next generalization step, which was 
made in [3], namely the introduction of discrete Toda type systems on arbitrary graphs. 

Definition 1. Let Q be a graph, with the set of vertices V{Q) and the set of edges E(Q). 
A discrete Toda type system on Q for a function x : V(Q) — > C reads: 



There is one equation for every vertex vo € V{Q); the summation is extended over 
star(uo), the set of vertices of Q connected to vq by an edge (see Figure [2]); we write 
xq = x(vq) and x = x(v) and often suppress the notational difference between the 
vertices v of the graph and the fields x = x(v) assigned to them. Often, the function 
cj) = 4>(xo,x;a) is supposed to additionally depend on some parameters a : E{Q) — > C, 
assigned to the edges of Q. 

The notion of integrability of discrete Toda type systems is not well established 
yet. We discuss here a definition based on the notion of the discrete zero curvature 
representation which works under an additional assumption about the graph Q. Namely, 
it has to come from a strongly regular polytopal cell decomposition of an oriented surface. 

We consider, in somewhat more detail, the dual graph (cell decomposition) Q* . Each 
t € E(Q) separates two faces of G, which in turn correspond to two vertices of Q* . A 
path between these two vertices is then declared the edge e* G E(Q*) dual to e. If one 
assigns a direction to an edge e £ E(G), then it will be assumed that the dual edge 
e* € E(Q*) is also directed, in a way consistent with the orientation of the underlying 
surface, namely so that the pair (e, e*) is positively oriented at its crossing point. This 
orientation convention implies that e** = — e. Finally, the faces of Q* are in a one-to-one 
correspondence with the vertices of Q: if xq € and xi,...,x n € V{Q) are its 

neighbors connected with xo by the edges ei = (xq, x\ ),..., t n = (xo,x n ) £ E(Q), then 
the face of Q* dual to xq is bounded by the dual edges e* = (yi, 2/2), • • • , ?n = (Un,yi)] 
see Figure El 





(26) 



star(?;o) 



7 




Figure 2: Star of a vertex x i n the Figure 3: Face of Q* dual to a vertex 

graph Q. xq of Q. 

We will say that a discrete Toda type system on Q possesses a discrete zero curvature 
representation if there is a collection of matrices L(e*; A) G G[X] from some loop group 
G[X], associated to directed edges e* G E(Q*) of the rfua/ graph Q*, such that: 

• the matrix L(e*; A) = L(xo,£,a; A) depends on the fields xq and a; at the vertices 
of the edge e = (xq,x) G E{Q), dual to the edge e* G E(Q*), as well as on the 
parameter a = a(e); 

• for any directed edge e* = (2/1,2/2)) if — e = (2/2,2/1)) then 

L(-e,A) = (L(e,A))- 1 ; (27) 

• for any closed path of directed edges 

*1 = (2/1)2/2), ^2 = (2/2,2/3), Cn = (2/n,yi), 

we have 

L(e;,A).-L(e;,A)^,A) = l. (28) 

The matrix L(e*; A) is interpreted as a transition matrix along the edge e* G E(Q*), that 
is, a transition across the edge e G E(Q). 

Under conditions (i27|) . ([28]) one can define a wave function ^ : V[Q*) — > G[X] on 
the vertices of the dual graph Q* , by the following requirement: for any directed edge 
e* = (2/1,2/2), the values of the wave functions at its ends must be connected via 

tf(y 2) A) = L(e*,A)v[/( m ,A). (29) 

For an arbitrary graph, the analytical consequences of the zero curvature represen- 
tation for a given collection of equations are not clear. However, in the case of regular 
graphs, like the square lattice or the regular triangular lattice, such a representation 
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may be used to determine conserved quantities for suitably defined Cauchy problems, 
as well as to apply powerful analytical methods for finding concrete solutions. 

It was shown in [3] that discrete Toda type systems with the following functions (f> 
are integrable in the above sense: 

(j)(x ,x;a) = , (30) 

X — Xq 

4>(xo,x;a) = acoth(x — xo), (31) 

<f>(x ,x;a) = log- — X ° + a , (32) 
x — xq — a 

, . . sinh(x-x + a) 

(p{x ,x;a) = log— r. (33) 

smn(x — xq — a) 

See [3] for details about the admissible assignments of edge parameters a. Actually the 
discrete Toda system with the functions (|3ip is reduced to the one with the functions 
(|30p via the change of variables x h- ► exp(2x) and therefore actually does not need to be 
considered separately. 



5 Quad-graphs and quad- equations 

Although one can consider 2D integrable systems on very different kinds of graphs on 
surfaces, there is one kind — quad-graphs — supporting the most fundamental integrable 
systems. 

Definition 2. A quad-graph T> is a strongly regular polytopal cell decomposition of a 
surface with all quadrilateral faces. 

Quad-graphs are privileged because from an arbitrary strongly regular polytopal 
cell decomposition Q one can produce a certain quad-graph V, called the double of Q. 
The double T> is a quad-graph, constructed from Q and its dual Q* as follows. The 
set of vertices of the double T> is V(T>) = V(G) U V(Q*). Each pair of dual edges, say 
e = (xo,xi) E E(Q) and e* = (2/1,2/2) £ E{Q*), defines a quadrilateral (xo, yi, x±, 2/2)- 
These quadrilaterals constitute the faces of a cell decomposition (quad-graph) V. Thus, 
a star of a vertex xq £ V(Q) generates a flower of adjacent quadrilaterals from F(T>) 
around xo; see Figure |U Let us stress that edges of T> belong neither to E{Q) nor to 

E{g*). 

Quad-graphs V coming as doubles are bipartite: the set V{T>) may be decomposed 
into two complementary halves, V(T>) = V{G) U V(Q*) ("black" and "white" vertices), 
such that the ends of each edge from E(T>) are of different colors. Equivalently, any 
closed loop consisting of edges of D has an even length. 

The construction of the double can be reversed. Start with a bipartite quad-graph 
T>. For instance, any quad-graph embedded in a plane or in an open disc is automat- 
ically bipartite. Any bipartite quad-graph produces two dual polytopal (in general, no 
more quadrilateral) cell decompositions Q and Q* , with V{Q) containing all the "black" 
vertices of T> and V{Q*) containing all the "white" ones, and edges of Q (resp. of Q*) 
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connecting "black" (resp. "white") vertices along the diagonals of each face of T>. The 
decomposition of V(T>) into V(Q) and V(Q*) is unique, up to interchanging the roles of 
G and G*. 

A privileged role played by the quad-graphs is reflected in the privileged role played 
in the theory of discrete integrable systems by the so called quad- equations supported 
by quad- graphs. 

Definition 3. For a given bipartite quad-graph T>, the system of quad- equations for a 
function x : V(V) — > C consists of equations of the type 

Q(x ,yi,x 1 ,y 2 ) = 0; (34) 

see Figure El There is one equation for every face (xo, yi, x%, 1/2) of T>. The function 
Q is supposed to be multi-affine, i.e., a polynomial of degree < 1 in each argument, so 
that equation (j3~4"|) is uniquely solvable for any of its arguments. Often, it is supposed 
that the function Q = Q(xo,yi,xi,y 2 ;a, (3) additionally depends on some parameters 
usually assigned to the edges of the quadrilaterals, a : E(D) — ► C, so that the opposite 
edges carry equal parameters: a = a(xo,yi) = a{y2,x\) and (3 = a(xg,y2) = a(yi,xi). 

There exists a fundamental and surprisingly simple notion of 3D consistency of quad- 
equations which can be put into the basis of the integrability theory, which has been done 
in [TO] and [13] . The property of 3D consistency allows one, in particular, to derive in 
an algorithmic way such basic integrability attributes as discrete zero curvature repre- 
sentations and Backhand transformations for quad-equations. Moreover, this property 
has been put [1] into the basis of a classification of integrable quad-equations which 
provided a finite list of such equations known nowadays as the "ABS list" . 

6 From quad-equations to discrete Toda type systems 

The geometric relation of a given surface graph Q to its double T>, described in Section 
[5j leads to a relation of discrete Toda type systems on Q to quad-equations on T>. The 
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2/2 



2/2 




latter relation is based on a deep and somewhat mysterious property of quad-equations 
which was discovered in several examples in [ID], was established for all equations of the 
ABS list in [3J, and was proved for all quad-equations with multi-affine functions Q by 
V. Adler, see Exercise 6.16 in [11] . 

Definition 4. A quad-equation (|34p possesses a three-leg form centered at the vertex 
xq if it is equivalent to the equation 

~4>(xo,yi) ~ ^0,2/2) = (p(x ,xi) (35) 

with some functions ip, (p. The terms on the left-hand side correspond to the "short" 
legs (xo, 2/1), (xq, 2/2) & E(T>), while the right-hand side corresponds to the "long" leg 

{x , Xl )eE(g). 

Summation of quad-graph equations for the flower of quadrilaterals adjacent to the 
"black" vertex xq G V{Q) (see Figure H]) immediately leads, due to the telescoping effect, 
to the following statement. 

Theorem 1. a) Suppose that equation (i34"|) on a bipartite quad-graph V possesses a 
three-leg form. Then the restriction of any solution f : V(T>) — * C to the "black" 
vertices V(G) satisfies the discrete Toda type equations, 

y~] 4>(x ,x k ) = 0. (36) 

star(xo) 

b) Conversely, given a solution f : V{Q) — > C of the Toda type equations (j36|) on 
a simply connected surface graph Q, there exists a one-parameter family of extensions 
f : V(T>) — > C satisfying equation (|34[) on the double T>. Such an extension is uniquely 
determined by the value at one arbitrary vertex ofV(Q*). 

It was shown in [ID] that symmetric discrete Toda type systems mentioned at the 
end of Section [3J come, through this construction, from the following integrable quad- 
equations: the systems with legs (|3P|) and ([32]) come from the 5 = and 5 = 1 cases, 
respectively, of the so called Ql equation of the ABS list, which reads 

a(x yi + xxy 2 ) - /3(x y2 + £12/1) - (a - (3)(x xi + yiy 2 ) + 5a(3(a - 0) = 0, 
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while the system with legs (|33p comes from the 5 = case of the so called Q3 equation 
of the ABS list, which reads 



sinh(a)(x yi + xiy 2 ) - sinh(/3)(x y2 + xiyi) - sinh(a - P)(x x 1 + yiy 2 ) = 0. 

7 Triangular lattice 

In Section [6] we established, for symmetric discrete Toda systems on an arbitrary planar 
graph Q, a relation to integrable quad-equations on the double T>. For non-symmetric 
discrete relativistic Toda type systems, such a relation remained unknown until recently, 
and it constitutes the main new result of the present paper. 

The non-symmetric discrete relativistic Toda type systems live on the regular tri- 
angular lattice T and cannot be directly generalized to arbitrary graphs. Therefore, 
we introduce now the specific notation tailored for the regular triangular lattice. The 
double of T is the quad-graph /C known as the dual kagome lattice (drawn on Figure [7] 
in dashed lines). The latter graph has vertices of two kinds, black vertices of valence 6 
and white vertices of valence 3, and edges of three types, all edges of each type being 
parallel. The quadrilateral faces of the dual kagome lattice are of three different types. 
We will denote them by type I, II, and III, according to Figure [H 



The dual kagome lattice can be realized as a quad-surface in Z 3 , so that the three 
types of quadrilaterals are realized as elementary squares of Z 3 parallel to the three 
coordinate planes (this is easy to see directly but follows also from the general theory 
of quasi-crystallic quad-graphs in [9]). In this realization, the black vertices of IC, that 



Xk-l 




Xk-l 



Xk+l 



Figure 7: Fields and wave functions on the triangular lattice 
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is, the vertices of T, are the points (^1,^2^3) £ Z 3 lying in the plane i\ + %i + 23 
while the white vertices of K are the points of Z 3 lying in the planes i\ + ^2 + 23 
(the vertices ^) and ii + 22 + 23 = —1 (the vertices <£). See Figure El 



= 0, 

= 1 




8 Discrete relativistic Toda type system from quad- equations 

on K, 

We now formulate the main result of the present paper. 

Theorem 2. Each discrete relativistic Toda type system is a restriction to the triangular 
lattice T of a certain 3D consistent system of quad- equations on the dual kagome lattice 
K, considered as a quad-surface in Z 3 . 

Proof of this theorem is obtained by a direct case-by-case construction of the corre- 
sponding systems of quad-equations (see, however, about the unifying "master system" 
in Section [9]). For the lack of space, these systems are given below not for all discrete 
relativistic Toda systems, but for four of them only, namely, for (fl~5l) . (fl6j) . (fT8l) . and 
(|20p . Details for other systems can be found in [12] . The systems are specified by giving 
the quad-equations explicitly for each type of quadrilateral faces separately in notation 
of Figure [9l One has to: a) find the three-leg forms, centered at x^, of quad-equations 
for all six quadrilaterals around and then check that adding these three- leg forms 
results in the corresponding discrete Toda equation, and b) check the 3D consistency of 
the quad-equations. All this is a matter of direct computations which are easy enough 
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to perform by hands but are better delegated to a symbolic manipulator like Mapl 
Mathematica. 

System (1151) . 3D consistent system of quad-equations: 

ahX(XY + UV) - ah(a - h)X 2 XU - ((a - h)X 2 + h) XV + a\ 2 YV = 0, 
aX(XY + UV) -XV + X 2 YV = 0, 

a(a - h)X(XY + E/V) + a/i(a - h)XU + ((a - h)X 2 + h) XV - aYV = 0. 
The three- leg forms of these equations (centered at Xk) read: 

( ae x k ~x k _ {a _ h) \ . g fo^ + A^ = ^ 

^ ' - (a - h) Xe x k e x k - aX® k+ i A' 

1 e Xk - aA$ fc +i Ae Xfc + a^k+i _ , 

1 e Xfc - (a - /i) A$ fc+ i * fc - aAe Xfc A 



ae x k ~ x k _ ( a _ h ) Xe x k + fr$ fe+1 ^ k h' 



(l + a 2 e Xk ~ Xk ~ l 



^ k -aXe x k ae x k + X$ k 
1 ae Xk + A$ fe $ fc - (a - /i) Ae Xfc 



l + a(a-/i)e x fc- x fc-i he Xk +X^ k 

Multiplying these equations leads to (fl~5|) . 

System (116ft . 3D consistent system of quad-equations: 

h(XY + UV) +YV-(1- hX) XV + h 2 XU = 0, 
a(IF + UV) +YV-(1- aX) XV + a 2 XU = 0, 
(/» - a) (AY + E7V) + (1 - aA)YF - (1 - hA)XY+ 
+ /i 2 (l - aA)XC7 - a 2 (l - hX)YU = 0. 



1. (I 
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\ U = e x k 



(a) North-western quadrilateral of type 
III 



Y = $* 



• Y = e Xk 



o 



/ U = $ fc+1 
'x = e x k 
(b) Northern quadrilateral of type I 



Y = $ 



fc+i 



h» JJ = e x k 



(c) Western quadrilateral of type II 



y = e x * 



^° x = * fe+1 

(d) Eastern quadrilateral of type II 




(e) Southern quadrilateral of type I (f) South-eastern quadrilateral of type 

III 

Figure 9: Notation for single quadrilaterals of the dual kagome lattice around the vertex 
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Three-leg forms of these equations, centered at x k : 



e x k -x k + h^k+l _ - h\)* k = Q 
e Xk + he x k 

_ ae x k+1 -x k _ £a>h + (1 - aXm+i = Q 

e X k e X k _|_ Q/\]/ fc+1 

(a - /Qeg+i ~^ (1 - aA)fr fc+1 | (1 - /iA)d> fc+1 _ Q 
1 - hae^ k + 1 ~ Xk e x k + a^ k+1 e x k + h$ k+1 
e x k -x k (l~hX)e Xk ^ 

*fc 3> fc + ae Xk 
(a - h)e x k~ x k-i | (l-qA)e x fc (1 - hX)e Xk ^ 
I _ hae Xk ~ Xk - 1 ® k + ae x fc + he Xk 

Adding these equations leads to ([TO . 

System (I18p . 3D consistent system of quad-equations: 

hX(XY + UV) - h\ 2 XU - (A 2 +h)XV + X 2 YV = 0, 
aX(XY + UV) - aX 2 XU - (A 2 + a) XV + X 2 YV = 0, 
{h - a)X{XY + UV) - h{X 2 + a)XU - (A 2 + h)XV+ 
+ a(A 2 + h)YU + (A 2 + a)YV = 0. 

Three-leg forms of these equations, centered at x k : 

A**. + he Xk e Xk h 



e x k — x k _ 



- Xe Xk e Xk - X§ k+1 X 
e x k _ \§ k+1 \ e x k _|_ a ^ k+ i 



l + ae Xk +i~ Xk e Xk e Xk - X^ k+1 

1 + ae% k + 1 ~ Xk e Xk - X^ k+1 Xe Xk + h$ k+1 
1 + heS k +i ~ Xk ' Xe Xk + aV k+1 ' e Xk - X$ k+1 

1 e Xk - X& k+1 - Xe Xk _ X 

e x k -x k _ 1 ' \ e x k + h<$> k+ i ¥~ k ~~ h ' 

r -r A $fc - Xe Xk 1 

1 _|_ Q. e x n X n -i \ _ » . jc _ _ 

/ Vk- Xe Xk \$ k + ae Xk X 
1 + he Xk ~ Xk - 1 A$ fc + ae Xfc - Xe Xk 



X, 

= 1, 



Multiplying these equations leads to (fTSj) . 



1. 
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System (1201) . 3D consistent system of quad-equations: 

h(XY -XU -YV + UV) - (1 - h\)(X — Y) — hX{U -V)- hX 2 = 0, (I) 

a{XY -XU -YV + UV) - (1 - aX)(X — Y) — aX(U -V)- a\ 2 = 0, (II) 
(h - a)(XY + UV) - h(l - 2aX)(XU + YV) + a(l - 2hX)(XV + YU)- 

- (1 - (a + h)X) (X — Y) — (h — a)X(U -V) - (h - a)X 2 = 0. (Ill) 

Three-leg forms of these equations, centered at x k : 

g t -*»)- 1 + *<*'--V> T ? = < N > 

(a* - $ fe +i - A) = 1, (E) 



l + a(x k+1 -x k ) ' " '"" ' x fc -^ fc+ i-A 
1 + a( x k+1 - x k ) x k - ^ fc+ i - A 1 - h(x k - + A) 
1 + h(x k+1 -x k ) 1 - a(x k - + A) x fc - $ k+1 - A 



(SE) 



Xfc " <i>fc+1 - A (x k -* k + X) = -l (S) 



x k - x k 1 - /i(x fc - ^>fc+i + A) " " /i ' 

1 + a (x k - x fc _i)) — ■ ~ - = 1, W 

x k -V k + X i + a (x k - $ k - X) 

1 + h(x k - x k -i) _ 1 + a(x fc - <S> fc - A) _ x fc - + A = ^ (nmv) 
1 + a(x k - x k -i) x k - + A 1 + h(x k - - A) 

Multiplying these equations leads to ([201) . 



9 The master system 

It turns out that all the 3D consistent systems of quad-equations leading to non- 
symmetric discrete relativistic Toda systems (those given above and those omitted for 
the space reasons), as well as the systems Ql and Q3,5 = o leading to the symmetric dis- 
crete relativistic Toda systems are particular or limiting cases of one multi-parametric 
system. Thus this latter system can be seen as the master one behind the whole the- 
ory of the relativistic Toda systems of the type ([Till (with discrete time) and ([I]) (with 
continuous time). 

Master system of quad-equations: 

{5 - (3-/)\(XY + UV) + (35{X 2 - >y)XU + {f3X 2 - 5)XV+ 

+ j(PX 2 -5)YU + (X 2 -j)YV = 0, (I) 

(r? - 7 e)A(IF + UV) + (A 2 - j)XU + (eA 2 - t])XV+ 

+ ^{eX 2 -rj)YU + eT](X 2 -j)YV = 0, (II) 

(/3rj - 5e)\{XY + UV) - (35{eX 2 - t])XU + erj(f3X 2 - 5)XV+ 

+ (f3X 2 - 5)YU - (eA 2 - rj)YV = 0. (Ill) 
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Three-leg forms of these equations, centered at x k : 



e x k + p e x k § e x k + A^ fc \ e x k - 7 j> fc+1 _ i 
7 e^fc+i 4. ^Xk e Xk _ A$ fc +i eAe Xfc + ^ k +i 



e Xk+i + ee x fc XeXk _ 7 $ fc+1 r^e^fe + A^ fc+ i 
p e Xk+i _ ge^fe ^gXfc + \\& k+1 \ e Xk + <5<j> fc+1 
§ e xk+i _ ^gXfc gAe^fe + qt k+1 e x k + /3A$fc +i 
7e Xfc + 5e~ k e Xk + /3A$fc +1 Ae Xfc - _ 
e x fc + (5e~ k Xe Xk + 5& k +i je x k - \q> k 
e x k + ee^fc- 1 je Xk - A^ fc Ae Xfc + r/^ k | 
7e x fe + ^Xfc-i Ae x ^ - e ^fc + e A$ fc " 
(Je^fe _ rj e Xk-i e Xk + e A$ fc /3Ae Xfc + ^ j 
^e^fc _ ee Xk-i \e x k + 77$*. J e x fc + A* k 



(E) 
(SE) 
(S) 
(W) 
(NW) 



Multiplying these three-leg forms leads to the following general equation of the discrete 
relativistic Toda type: 

e Xk~Xk _|_ p ^ e Xk~Xk_i_§ 



<y e X k -X k+ § e X k -X k+ p 

-yg^-fc+l X k _|_ „ ^X k IH 1 f 



_|_ g ^ys x k 1 _|_ ^ 

f3e~ k + 1 ~ Xk - e §e Xk ~ x k-i 



Se%k+1 x k-jj p e x k -x k -i 



1. (37) 



This equation has five parameters /?, 7, <5, e, 77. Actually there are only four, because of 
homogeneity: if 7 ^ 0, we can set 7 = 1 by replacing 5, rj through 5/7, 77/7, respectively. 
Moreover, we could eliminate two further parameters by shifts of the form x k {t) — > 
x k (t) + Ak + which keep the form of the equation invariant. 

It is not difficult to find out the special values of parameters which lead to all the 
discrete relativistic Toda type equations listed in Section In particular, the value 
7 = 1 leads to the most general symmetric equation (|22p . with further degenerations to 
(|23p . (|24p . (|25p . The value 7 = is of the primary interest for the aims of the present 
paper, since it leads to 

e x k —x k _|_ p e x k — x k _i _|_ e p e x k +\ — X k _ g $ e x k — X k ^\ _ ^ 
e x k -x k+ p e x k +i-Xk +e fi e Xk+i-x k _ r j p e x k -x k _i_ e 

which happens to encapsulate all the non-symmetric equations. For instance: 

• System (fT5|) appears from ([38]) with the choice rj = 00, f3 = (h— a) /a and e = 1/a 2 . 
In the quad-equations it is convenient to set 5 = ah and to re-scale A qA. 
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System (|18j) appears from (|38p with the choice /3 = rj = — 1, 5 = h and e = 1/a. 

One gets from (1381) to the additive equation (fT6l) in two steps. On the first step, 
one starts with the parameter values (3 = (h — 0)/6, r] = 9/(a — 9), 5 = h9, 
e = l/(a9), which leads to a remarkable equation introduced in [18J: 



l + 0h- 1 (e Xk ~ Xk - 

l + 9h- 1 (e Xk ~~ k -l) 

_ l + 9ae x k+i~ x k \ + h(9 -a)e~ k+1 ~ Xk 1 + a(9 - h)e x k~ x k-i 

l + 9ae x k~ x k-i l + a {9-h)e~ k + x ~ Xh 1 + h{9 - a) e x k~ x k-i 



• (*) 



This equation interpolates between (fT5j) (corresponding to 9 = a) and (fTUj) (which 
corresponds to 9 = 0). In the corresponding quad-equations it is convenient to 
re-scale A 9 A. We remark that the last equation is a time discretization of 

/ g^fc+i x k e x k x k —\ 

x k = (1 + 9x k ) (1 + ax k+1 ) (l + ax k -i) 



l + 9ae x k+i~ x k i + ae Xk~x k -i 

e 2(x k+ i-x k ) e 2(x fe -x fc _i) \ 

+a(9 — a) a(0 — a) , 

l + 9ae x k+i~ x k i + ae x k -x k ^ 1 J 

which in turn interpolates between the continuous time equations ([2]) (for 9 = a) 
and (J3]) (for 9 = 0). On the second step, one performs in equation Q the limit 
9^0. In this limit one should also re-scale the auxiliary variables according to 
q> ^ A^, $ ~> $/A, and A 1 + 9X/2. 



To transform (j38|) to the rational equation f)20|) one makes the change of variables 

x ~-> KX, 

accompanied by the change of parameters 

[3 —i -|- S ~* — 1 + k5, e —1 + ne, t) ~~> — 1 + fcrj. 
Sending k — > 0, one arrives at 

Xk - x k + f3 x k - Xfc-i + e x fc+ i - Xfc - j3 + e x fc - x k -i - 5 + rj _ 
x k - x k + j3 x k+ i - x k + e -x k - 5 + T] x k - x k -i - f3 + e 

Equation ([20]) corresponds to the choice 

/3 = rj = 0, 8 = -l/h, e = l/a. 
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10 Zero curvature representations 



The construction of discrete Toda type systems on graphs from systems of quad-equations 
allows one to find, in an algorithmic way, discrete zero curvature representations for the 
former. Indeed, each quad-equation can be viewed as a Mobius transformation of the 
field at one white vertex of a quad into the field at the other white vertex, with the 
coefficients dependent on the fields at the both black vertices. The SL(2, C) matri- 
ces representing these Mobius transformations play then the role of transition matrices 
across the edges connecting the black vertices. The property (|28|) is satisfied automati- 
cally, by construction. 

Specializing this construction to the case of the regular triangular lattice (see Figure 
[7]), we denote by L k the transition matrix from ^f k to ^fe+i, and by Vk the transition 
matrix from ty k to ^> k . I n this notation, the discrete zero representation reads: 

L k V k = V k+1 L k , (39) 

both parts representing the transition from ^ k to *& k +i along two different paths. It is 
clear that L k is the product of two matrices, the first corresponding to the transition 
from ^ k to 3>fc+i across the edge [x k , x k ], and the second corresponding to the transition 
from ^fc+i to ^fc+i across the edge [x k , x k+ i], so that 

L k = L(x k , x k , Xk+r, A). (40) 

Similarly, V k can be represented as the product of two matrices, the first corresponding 
to the transition from ^ k to across the edge [x k , Xk-i\, and the second corresponding 
to the transition from & k to ^fr k across the edge [x k ,x k _i], so that 

V k = V(x k ,Xk-i,x k -i; A). (41) 

The matrices L k , V k for a given discrete relativistic Toda type equation can be computed 
in a straightforward way, as soon as the generating system of quad-equations is known. 
Theorem [2] provides us with the means for this goal. 

The resulting zero curvature representations possess an additional remarkable prop- 
erty. It is well known (see, e.g., [3], [J]) that the discrete relativistic Toda type equations 
possess a Lagrangian (variational) interpretation with a discrete Lagrange function 

£( x > S) = X] ( Al ( x fc ~ Sk) ~ A 2(x k+1 - x k ) - J^z{x k+ i - £ fc )) • ( 42 ) 
fcez 

Here Ai, A2, A3 are antiderivatives of the functions F, G, H in the general equation (|14|) . 
The canonically conjugate momenta and the Lagrangian form of equations of motion 
are given by 

dC(x, x) _ dC(x,x) 

Pk 7, 7. ) 

ox k dx k 

which specializes in our case to 

p k = F(x k - Xk) + H(x k+ i - x k ) (43) 
= F(x k - x k ) + H(x k - x k -i) - G(x k+ i - x k ) + G(x k - x k -i). (44) 
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Theorem 3. For all discrete relativistic Toda type systems, the transition matrix L k 
from equation [JIfy is local, when expressed in terms of canonically conjugate variables: 



L k = L(x k ,p k ;X). (45) 

Moreover, as a matter of fact, the matrix L k does not depend on the time discretiza- 
tion parameter h, so that the corresponding Lagrangian maps (x,p) \— » (x,p) belong to 
the same integrable hierarchies as their respective continuous time Hamiltonian coun- 
terparts. In other words, these Lagrangian maps serve as Backlund transformations for 
the respective Hamiltonian flows, the Backlund parameter being the time step h. 

Proof of this theorem is obtained again via direct computations on the case-by-case 
basis. For all cases where the local discrete zero curvature representation was known 
(those cases are listed in [18]) the results obtained from the system of quad-equations 
coincide with the previously available ones. Therefore, we illustrate the claims of the 
theorem with the case where the local discrete zero curvature representation was not 
known previously, namely with the rational systems (|20f) . (121j) . It is useful to keep in 
mind that these two systems come as two elementary flows (a positive and a negative 
ones) of the same hierarchy [18] . □ 

Equation ( 1201) . The Lagrangian form reads: 

= (x k -x k ). ) + ^ k+1 - Xk \ (46) 
1 + a{x k+1 - x k ) 

( ~ _ \ 1 + a{x k - x k -i) 1 + h(x k - x k -i) 

— ( x k x k) ' 1 / \ ' 1 i / ~ \ ' l 4 V 

1 + a{x k+ i - x k ) 1 + a{x k - x k -i) 

The transition matrices of the zero curvature representation of this map read: 

f-X + x k X 2 + XaePk - (x k - aePk)x k - ePk\ 
Lk= ( 1 -X-x k + aeP* J (48) 



and 




A 2 + X(x k - x fc _i + aePk-i) - (x k _i - aePk-i)x k 
-X - x k -\ + aePk- 1 



(49) 



Note that in the limit h — * one obtains a zero curvature representation ([7]) for the 
Hamiltonian form of equation ([8|) with the same matrix L k as in (|48p and with the 
matrix 

_ f-X + x k A 2 + X(x k - x fc _i + aePk-i) - (x k -i - aeP*>-i)x k \ . , 
Mk ~ V 1 -A - + aePk~i J ■ {bU) 

To the best of our knowledge, this result (and even its non-relativistic particular case 
a = 0) was previously unknown. 
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Equation (1211) . The Lagrangian form reads: 

^ = l + a(a + h)h-Hx k -x k ) ' (1 + ( ° + h){Xk ~ ~ Xk ~ l]) (51) 
x k~x k l + a(x k+1 -x k ) 

■ (1 + (a + h)(x k+1 - Xk)). 



1 + a (a + h)h~ 1 (x k — x k ) 1 + a(x k — Xk-l) 

(52) 

This Lagrangian map admits a discrete zero curvature representation 

W k+1 L k = L k W k (53) 
with the same transition matrix L k as in (|48p and with 

ww * 



(1 - 2aX) i 1 + (a + /i)(x fc - - aePfc 

A + x fc -ae^ fc A 2 + X(x k — x k -\ — aeP k ) — (x k — aeP k )x k -\\ , . 

1 A - / 

Again, in the limit h — > one obtains a zero curvature representation ([7|) for the Hamil- 
tonian form of equation with the same matrix L k as in (|48l) and with the matrix 



M fc = ? s - 

(1 - 2aA) [1 + a{x k - x fc _i - aePfc)J 

(55) 



1 A - Zfc-i 

Also these results seem to be previously unknown, even in the continuous time case 



11 Conclusions 

In this paper, we clarified the origin of all non-symmetric discrete equations of the rela- 
tivistic Toda type from 3D consistent systems of quad-equations. Unlike the symmetric 
case, the three coordinate planes carry here three different quad-equations, so that a 
more general understanding of the 3D consistency than usual is required. Note that this 
more general concept has already been discussed and laid into the basis of a classifica- 
tion procedure in [5], see also examples discussed in [8]. However, the classification only 
has been performed for the so called systems of type Q in [5] (all edge biquadratics non- 
degenerate, see the original paper for details). Examples which arose in the present work 
demonstrate that also the systems of type H (with some of the edge biquadratics being 
degenerate) are of a considerable importance, which calls for a complete classification 
of this case, as well. We plan to turn to this tedious task in our future work. 
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